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Abstract 

In this paper we consider 2+1-dimensional gravity coupled to N point- 
particles. We introduce a gauge in which the z- and z-components of the 
dreibein field become holomorphic and anti-holomorphic respectively. As 
a result we can restrict ourselves to the complex plane. Next we show 
that solving the dreibein- field: e^{z) is equivalent to solving the Riemann- 
Hilbert problem for the group 50(2, 1). We give the explicit solution for 
2 particles in terms of hypergeometric functions. In the N-particle case we 
give a representation in terms of conformal field theory. The dreibeins are 
expressed as correlators of 2 free fermion fields and twistoperators at the 
position of the particles. 



^E-mail: welling@fys.ruu.nl 



It is well known that in 2+1 dimensions, Einstein gravity is a topological 
theory. One way of understanding this is by counting the degrees of freedom of 
the gravitational field. One finds 6 independent metric components, 3 first class 
constraints which generate 3 gauge transformations (coordinate transformations). 
These have to be fixed by gauge conditions, introducing 3 more constraints. So 
the gravitational field has no local degrees of freedom (there are no gravitational 
waves or gravitons). It is however possible to introduce some degrees of freedom 
by the introduction of point particles. In the case all these particles are static, 
the metric tensor ga/s was solved by Deser, Jackiw and 't Hooft 0. Also a 
geometrical approach was invented by 't Hooft to treat moving point particles 

0. Here it was proved that a configuration of moving point particles admits a 
Cauchy formulation within which no closed timelike curves (CTC) are generated. 
Spinning particles however do generate CTC's and this is the reason why we 
will not consider them in this paper. Quantum-mechanically one hopes that this 
problem will be resolved. 

To attack problems in 2+1 dimensional gravity there are 2 mainstreams. 
There are researchers who use Wittens Chern-Simons formulation 0] and there 
are researchers who stick to the ADM-formalism. In the CS-formulation, point 
particles are coupled to non-abelian gauge fields (that take values in the group- 
algebra iso{2, 1)) by introducing " Poincare-coordinates" q"'{x), that transform as 
poincare vectors under local gauge transformations We feel that the geomet- 
rical meaning of the variables in the Chern-Simons approach is less clear than in 
the ADM formulation. We will therefore stick to the ADM-formulation of gravity. 
Furthermore we will only consider the case of an open universe. 

Although the classical theory of gravity in 2+1 dimensions is relatively well 
understood, the first-quantized models are not all equivalent and it is not 
certain that a consistent quantum theory exists. To our opinion however it is still 
interesting to find out what produces these problems, and see if it is connected 
with the fundamentals of Einstein-gravity (covariance) or if it is just an artifact 
of the 2+1-dimensional model. It might be especially interesting to see if one can 
find a consistent S-matrix for scattering of quantum-particles. The scattering of 
2 quantized point particles was first studied by 't Hooft and later by Deser and 
Jackiw 10. In the CS-formulation scattering of particles was studied by Carlip 
and later by Koehler et al.p|. To our knowledge a second-quantized model for 

1. e. a scalar field coupled to gravity has not yet been investigated. In such a 
model one could study the ultraviolet behaviour of Einstein-gravity because 
one expects to have much better understanding of the non-perturbative theory 
than in 3+1 dimensions. 

In what follows we will restrict ourselves to the classical theory of gravity in 
2+1 dimensions, coupled to a set of moving point particles. The total mass of 

^In 2+1 dimensions gravity coupled to a scalar field is still non-renormalisable because the 
coupling constant is dimensionfuU. 
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the particles must be such that the universe does not close. This is necessary for 
our gauge choice to be valid for all times t. In this case we also avoid the problem 
of time, because we have a boudary term that will act as a Hamiltonian. So we 
use the asymptotic coordinate t as our time- variable. 

In this paper we start to review some well known aspects of general relativity 
in section 1. It is stressed that for an open universe we need a boundary term 
in the action which serves as a Hamiltonian for the reduced theory. The reduced 
theory must be obtained by putting the solutions for the metric g^j^y back into 
the action. In 2+1 dimensional gravity the gravitational field can in principle be 
completely integrated out because it has no dynamical degrees of freedom. In 
the next section we introduce a slicing condition in such a way that the area of 
the Cauchy surface is locally maximal. On that slice we then choose conformal 
coordinates. In this gauge the coordinates are not multivalued as in the "flat 
coordinate-system" (in the flat coordinate-system the metric is everywhere rjni, 
except at the particles positions) . Another pleasant feature of this gauge is that 
the tetrad field e^" defined by the relation g^^, — "e^r^ab, in complex coordinates, 
splits into a holomorphic and a anti-holomorphic sector. Actually all information 
is contained in the field e/, which is a holormophic function and can be studied on 
the complex plane. A residual gauge invariance are the conformal transformations 
and we can use them to bring 2 particles to definit positions (leaving infinity 
invariant). Another pleasant property of the coordinates is the fact that in the 
process of reducing the action, only the kinetic term of the matter variables 
together with the boundary term survives. As mentioned, the boundary term 
takes the role of the Hamiltonian for the matter fields. Next we solve in our 
coordinate system the gravitational field for a stationary, spinning particle with 
mass M, sitting at rest in the origin. Although this metric is pathological near the 
particle, the asymptotic behaviour of this solution is very useful! as asymptotic 
condition for the general solution. This makes sense because we expect that any 
configuration of particles with total mass M and total angular momentum J can 
be described by an effective center of mass particle with mass M and spin S=J. 
In section 3 we refrase the problem of finding solutions for the tetrads e/ into 
an old mathmatical problem: The Riemann-Hilbert problem. The Riemann- 
Hilbert problem deals with a vectorfield e^", in general living on a Riemann- 
surface with N singularities, that has prescribed monodromy properties. In our 
case we will work on the Riemann-sphere and the monodromy must take values 
in S0(2,l). Actually it proves more convenient use the spinor representation of 
SU(1,1) and find a spinor (^"{z) with the corresponding monodromy properties. 
The Riemann-Hilbert problem can be attacked in 2 different ways. The first 
method one can use seeks a d'th order linear differential equation with N regular 
singularities, (d is the dimension of the representation used; 3 for S0(2,l) and 
2 for SU(1,1)). Throughout the paper we stress the importance of the local 
behaviour of the solution near the singular points (particles) and at infinity. As 
an example of this method we solve the 2 particle case in terms of hypergeometric 
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functions. The N-particle case is more conveniently studied using a first order 
linear matrix differential equation where the dimension of the matrices is d. We 
use a representation due to Miwa, Sato and Jimbo to represent the solutions 
to this equation in terms of correlators of 2 free fermions with twistoperator 
insertions on the positions of the particles. We check that the local behaviour 
of this solution near the singular point indeed agrees with what one expects. As 
we move the singular point (for instance by evolution) the monodromy does not 
change. This implies that we can write down a set of Schlesinger deformation 
equations. They describe basically how the matrix differential equation and its 
fundamental matrix of (rationally independent) solutions change if we change the 
deformation parameters (i.e. the location of the particles). Next we check that 
we can reproduce the static solution of N particles at rest. It should be noted 
that this special example bears much resemblance with the Coulomb gas. Finally 
we derive some kind of Knizhnik-Zamalodchikov equation for the correlator of 
the twistfields (without the fermions). 
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1 The ADM- formalism 



We will work on space-times with the topology Ex i? where E is an open orientable 
2 dimensional surface. In the ADM-formalism one makes a 2+1-split of the metric 
in the following way: f\ 

= {-N^ + g,jN'N^)dt^ + 2Nidx'dt + Qijdx'dx^ (1) 

In this paper roman indices in the middel of the alphabet (i,j,k,...) will range 
over 1,2, roman indices in the beginning of the alphabet (a,b,c,...) are flat indices 
and run over 0,1,2 and greek indices are space-time indices and run over 0,1,2. 
The indices of 2 dimensional tensors are raised and lowered by gij. This object 
is also the dynamical variable of the theory. The lapse function and the shift 
functions Ni are Lagrange-multiplier fields and as such not dynamical. They 
determine the way the 2 dimensional surface is inbedded in the 3 dimensional 
space-time. The field canonically conjungate to gij is: 

^{K^^ - K\g^^) (2) 



TT 
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Here ^ is the square root of the determinant of the 2-dimensional metric, g'^^ is 
the inverse of g^j {gijg^'' = 5'), and Kij is the extrinsic curvature that measures 
how the normal to the surface S changes as we walk along that surface. It is 
determined by: 

Kr, = ^idtg.j~2D^,N,)) (3) 

Here Di is the covariant derivative with respect to the 2-dimensional surface 
defined by gij. The brackets around indices means that we symmetrize the ex- 
pression in those indices. 

The Lagrangian for gravity coupled to matter in terms of these variables is: 

(4) 
(5) 

(6) 

(7) 
(8) 

where means the intrinsic, 2-dimensional curvature scalar and x' = The 

' or 

equations of motion derived from this action are then the Einstein equations: 

G^, = 8nGT^, (9) 



L = 


Lg + Lm 




= 


j d?x Tf'^dtgij - 


- NH - N'Hi 


H = 




Y) - V7 ^'^R 


= 






Lm = 


-IQttG f dTj2 

•' n 


Tn„ \J —g^yXn x'^ 



^Most of the expressions in this section can be found in jl] 
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where 



rpfjLV 



The expression for the matter Lagrangian can also be written in a first order 
form: 

Lm = -IQttG J dr [p.e^'^x'" - ^{v^'PaPb + m')] (11) 

This form is more convenient for us as we will work with the "dreibein" e,". Of 
course the metric is easily recovered from the dreibein by the formula: 

9^.^ = e^^'ej'rjab (12) 

Pa Has the interpretation as the momentum of the particle in a fiat coordinate 
system. One immediately recognizes what a formidable task it is to solve these 
equations in the general case of N moving particles in arbitrary directions with 
arbitrary velocities. We will therefore attack the problem from a different per- 
spective. 

Finally we want to remark that in the case of an open universe Regge and Teit- 
elboim [|I^ found that one should add to the Lagrangian (^) a surface term that 
is needed to have a well defined variational derivative. This surface term is then 
identified with the total energy contained in the universe. In 2+1- dimensions 
this surface term is 0: 

It is now interesting to perform a reduction of the system to the physical 
degrees of freedom. One has to insert the solution of the gravitational fields into 
the total Lagrangian (^. In the coordinates chosen by us, this implies that Lg 
will vanish completely and the surface term becomes the Hamiltonian for the 
reduced system. So it is an "effective" Hamiltonian that describes how the point 
particles interact. The fact that there are no terms that survive in Lg is clear if 
one remembers that the gravitational field itself contains no degrees of freedom. 
We will return to this point in the discussion. 



2 Gauge fixing and asymptotic conditions 

The problem of deriving the gravitational fields around moving point particles is 
already solved in the fiat coordinate system |Q in a geometrical approach. The 
fiat coordinates tt" are however multivalued. One can see this already in the 
simplest configuration: a particle sitting at rest in the origin. The space-time 
is constructed by cutting a wedge out of space-time. Consider now a point 
somewhere in the plane (see fig.p. If an observer approaches the point B"^ along 
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y-axis 




Figure 1: Space-time around a particle sitting at rest in the origin 

path Pi he will find coordinates: 

^„^/^cos^\ 

However, an observer moving along path /?2 would give point B°' the coordinates: 

/ i?cos(</? + a) \ 

\ Rsmi^ + a) ) ^^^^ 

One observes that it is really the unusual ranges (or matching conditions) of 
the coordinates that make them multivalued. The main aim of this paper is to 
introduce new, single- valued, coordinates {z,z), that transform the fiat metric 
Tjab into a nontrivial metric Qa/s according to the well known transformation law: 

g,j,u = d^.u'^dr.u^-qab (16) 

Part of the gauge fixing-procedure is choosing a Cauchy-surface. The Cauchy 
surfaces we are looking for, are conveniently characterized by the property that 
they locally maximize the area of the surface. Since we know that locally there 
must always be a surface that has maximal area, we know that a solution to this 
coordinate condition always exists (i.e. the gauge is accessible). This is believed 
to be the simplest gauge condition that provides for single valued coordinates 
{z, z, t), even though the elaboration of this condition is fairly complicated as will 
become clear. The coordinates u°'{z, z; t) arc considered as embedding conditions 
for each time t. They are flat coordinates but unfortunately multivalued functions 
on these Cauchy surfaces. So in order to descibe them we need cuts in the surface, 
called strings in previous literature. The transformation of u"" over a string is in 
general a Poincare transformation. We can also define the so called dreibein or 
tetrad fields on the surface E: 

e; = (17) 
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They are also multivalued but transform by homogeneous Lorenztransformations 
if we move over a string. In what follows we will often work in complex coordi- 
nates: 

z = (x + iy) 

z = (x — iy) (18) 

It is easy to write down an action that would produce the slicing condition that 
the area of the Cauchy surface S is locally maximal. We can view the coordinates 
{z,z) as coordinates of a Euclidean worldsheet living in a fiat 3 dimensional 
target-space described by the coordinates m*^. The embedding conditions for the 
worldsheet to have maximal area are then precisely given by the Polyakov action: 

s 

Sp = j <fz Vl9''d,u''dju%, (19) 

where gij is the 2 dimensional metric on the worldsheet, 7 =det gij and rjab is the 
3 dimensional metric in the targetspace. We can also view the g^j as the induced 
metric on the Cauchy surface derived by projecting r]ab on S. If we furthermore 
choose these coordinates on the worldsheet to be conformally flat: 

g^J = (20) 

then this action reduces to: 

Sp = J (fz du''du\b (21) 

where we used the complex coordinates {z, z)) and real coordinates u"". Further- 
more we used the notation: ^ = d ^ = 9. It is well known that the equations 
that follow from the string action are: 

du^du^r^ab = (22) 
du''dv!'7]ab = (23) 
ddu" = (24) 

This implies that u"" is a harmonic function: u"'{z, z; t) = u°'{z; t)+u°'{z; t). These 
conditions imply for the dreibein field ([T7|) : 

(25) 

(26) 
(27) 
(28) 
(29) 





= C\t) + dt j dz e: + dt j dz 


e/ = du" 


= e^{z) (holomorphic) 


e/ = Bu" 


= 6/(2) (anti- holomorphic) 




= (null-vector) 




= (null-vector) 



"'in string theory one could say that the string follows a trajectory so that the worldsheet 
has maximal area. Because we consider a Euclidean "worldsheet" we can not speak about 
trajectories. 
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where C"(t) is some vector only depending on t. Note that we still have a remain- 
ing conformal freedom: z w{z). This will be used later to transform 2 particles 
to definit positions (leaving infinity invariant). Using equations (pOD,(p!2D,(|3|) and 
(0) one can calculate what the implications are for the fields gij, vr*-^. 

gij = eHij (30) 
tt^ttI = (31) 

The advantage of this gauge choice is now evident. First of all, as we will see 
later, we will only have to solve the holomorphic, nuU-vectorfields e/, living on 
the complex plane, in order to find the gravitational field gap- The fact that the 
dreibein splits into a holomorphic and a anti-holomorphic part is a peculiar fact of 
2+1 dimensional gravity. From a mathematical viewpoint it is a convenient gauge 
condition because we have now the machinery of complex calculus on Riemann 
surfaces at our disposal. The conformal freedom also suggests that there might 
be a connection with conformal field theory. As we will show later, the dreibein 
is formally solved as a two-point function of free conformal fermion fields with 
some twistoperator insertions. Another advantage is that in the reduction scheme 
explained in section |I], the reduced "effective" Lagrangian looks very simple. Since 
TT is the canonicacally conjugate momentum of 26*^ and we have formulas (^) and 
(|3TD we find that in the Lagrangian (^) all kinetic terms of the gravitational field 
vanish: This is easily seen by splitting the kinetic term as follows: 

ix'^dtg^, = in^Y'dtgJj + nOte'^ (32) 

where vr = 7r| and a traceless tensor P^^ is defined by 

= P^J - \Pl5^J (33) 

So because gj^ = and tt = we see that this kinetic term vanishes. Furthermore 
the term: NH + NiH^, which also contains matter contributions, must vanish 
because these are precisely the constraints that we solve. The reduced La- 
grangian is thus given by the kinetic term of the matter fields and the surface 
term (pTSj) which now only depends on the matter phase space variables because 
we (formally) solved the gravitational fields in terms of the matter variables. The 
surface term thus serves as a hamiltonian for the matter degrees of freedom. We 
will see later that it represents actually the total deficit angle of the universe (see 
[pnp. As an example, we will work out the solution found in for a massive, 
spinning particle at the origin, in our gauge. The solution for the metric will be- 
come degenerate near the particle. We are however not interested in that region 

^More precisely, the set of first class constraints {H, iJ'} becomes after imposing the gauge 
conditions (^0|) and (^) a twice as large set of second class constraints which can be set equal 
to strongly on the constraint surface. 
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since it is exactly the place where closed time-like curves may occur. At infinity 
this solution is however very interesting as we expect that any configuration of 
particles can be described asymptotically by an effective center of mass parti- 
cle with mass M equal to the total energy contained in the universe and spin S 
equal to the total angular momentum of the universe. The solution therefore will 
provide the asymptotic behaviour of the gravitational fields. 

If we use fiat coordinates u"" the matching conditions for this particle (sitting 
at rest in the origin) will be 0: 



pa 
JX u 




(34) 
(35) 

(36) 



Here u"" is the value of u"" after an anti-clockwise rotation about the particle. 
Where we have defined: a = 4GM and A = AGS. For a general choice of origin 
and a moving particle the identification will be an element of the Poincare group: 



(37) 



Here B\ is a boost matrix with arbitrary rapidity (r/) and in an arbitrary di- 
rection, and q"' is now a general translation (not necessarily of the form (0)). The 
line-element for the massive spinning particle sitting in the origin (not moving), 
found in p|, is given by: 



-{dt + AdOf + ^(rfr^ + r'^dd 



(3^ 



This is clearly not in the gauge we have chosen and to transform this into our 
gauge we have to write it as: 



ds" 



-dt^ - 2Adtd9{p, u) + e^{dp^ + p^du"^ 



(39) 



where the difference is in the A^dO'^ term. We took here also Qoo = — 1 because the 
system is stationary and goo = dtu"'dtu^r]ab- The map from the {t, r, 6')-coordinate 
system to the {t, p, u;)-coordinate system is provided by: 



In 



u 

P_ 
Po 



e 



1 — a 



arccosh(- 



A 



(40) 
(41) 
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If we choose tq^ " = A and po = p{'>^o) = 1 and use the complex coordinates: 
p = zz, u = Yi lii(f ) we find for the fine element in our gauge: 

ds"^ = -dt^ + —dtdz - —dzdt + —[^^^ 2 + A^{zz)"'^]dzdz (42) 

Z Z 4:ZZ l\ 

So this will be the asymptotic behaviour of any metric describing N moving point- 
particles with total energy ^ and total angular momentum From here one 
finds that the coordinates are given by: 

A A 

M° = t + — In z In z (43) 

u = ^^^^^ + 4—] (44) 

v^(l-a)^ ^1-"^ ^ ' 

u = r[^'"" + ^] (45) 

v^(l-a)^ 2;^-" 

We see that both the vP and the u and u are multivalued. Note that the log 

terms in vP are needed to generate the time jump of 871GS. The dreibein fields 
are given by: 

aoM" = eS = (1,0,0) (46) 

A 1 /12 a-2 

5,M" = e/ = ( — V ) (47) 

a,?/'^ = e/ = (-_,-ii^,_^) (48) 



one can easily check that these dreibein fields fullfill our gauge conditions (^ ,27,28,29) 
These coordinates m° and the 3 dreibeins are only determined up to a global 
Lorentztransformation, since A^^ti* would produce the exact same line element 
(^21). Actually we can choose our global Lorentzframe by the condition that the 
dreibeins have this asymptotic form. This immediately implies that the constant 
function C"(t) in ( P^D is determined to be 6q. In general however we will allow 
for more general global Lorentzframes (i.e. frames in which the center of mass 
moves). It will be this this global frame that determines the function C"(t). 



3 the Riemann-Hilbert problem 



The major advantage of our gauge choice is that the dreibein field splits into a 
holomorphic and a anti-holomorphic piece (and e^"'). In effect we only have to 
solve the e/-field because the e/ is then given by: 



n(e. 



F 




(49) 
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Then, in the Lorentzframe where the effective center of mass is not moving, the 
eo"-field is given by the the expression: 

Cg'^ = 5^ + dtJdz e/ + dt J dz e/ (50) 

So we only have to solve for the null-vector e/, living on the complex plane. 
We will not try to solve the Einstein-equations explicitly but we will use the 
monodromy properties of the dreibein. Now what happens if we parallel transport 
the dreibein ej^ over a loop 7„ around a point particle at position a„. First we 
notice that we can deform the loop 7„ as long as we don't pass through a particle. 
This means that the tranformation will only differ on different homotopy classes 
of loops. To be more precise, we are interested in the transformations from 
the first homotopy class of loops on a punctured Riemann surface, to the group 
S0(2,l): 

7Vi{C - {ai, ...aN, oo}; Zo) S0(2,l) (51) 

We would like to view the coordinates u"'{z) f\ as multivalued functions on the 
complex plane. We therefore need a cut in the surface where the functions ?/'*(z) 
transform. We attach to each particle a string over which the u"'{z) transform as 
a Poincare vector: 

u"" {BRB-^ ■ uY + l" (52) 

At this point we should warn the reader that the translation vector is not the 
same as the translation vector q"" in (P?]). This is due to the fact that we consider 
here only the holomorphic part of u"". For the anti-holomorphic part we have 
a similar transformation with a translation vector (l*)"". The correspondence is 
thus given by q"" = l"" + (/*)", which is a real vector. The transformation law 
(|52D implies that the e^" = du"" transform under the homogeneous part of the the 
Poincare group; the Lorentzgroup. So if we parallel transport the dreibein along 
a closed curve (loop) 7„ we find: 

e/ 'i" {B^RnB-' ■ e,r (53) 

As the particles evolve we don't allow them to hit a string, so one cannot avoid 
that the strings might get 'knotted' in a very complicated way. On the other 
hand, if we have an expression for we can reconstruct u°-{z) by integrating 
once: 

u^'iz) = I dz ej" 



^The vector u°'{z) is the holomorphic part of the vector u'^{z, z) — u°-{z) + (see section 

2). 
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Figure 2: Particles with strings attached moving on the complex plane. 



It is now easy to see that the multivaluedness of is indeed of the form (52). 
If we integrate around a singularity at we find: 

v!'{zo) = / dze;+ i dz e/ + Mk dz e/ (54) 
= (1 - Mk)Al + (f dz e/ (55) 

Jafe 

where = dz is the position of the particle in flat coordinates, and 
Mk = BkRkB^^ (the monodromy matrix). Because u{zo) = is the origin of 
the flat coordinate system, the u"" represent the shift of the origin and as such is 
equal to the vector /° of ([5^). If we integrate around 2 particles we find: 



u''{zo) = (1 - Mi)A1 + Mi(l - M2)Al + i dz e/ + f, i dz e^z) (56) 

From this expression, together with the antiholomorphic part, we can derive the 
relative angular momentum of the 2 particles. It is defined in the center of mass 
frame of the 2 particles by: 

L = [B-L. ■ {I + nf (57) 
This is precisely the same expression as in only if we demand: 



lim i dz e/ = (58) 

where we take the radius R of intgration to go to zero. To produce the multival- 
uedness of the e/ we expect, in the frame where the particle is boosted to rest, 
the following local behaviour near the location of a particle ai. 

{B-\,r'^'^ z'^Ht{z) (59) 

The Af (called local exponents) are real numbers and the H^{z) are holomorphic 
functions near Oj. The condition (p8|) now implies that: 

K > -1 (60) 
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In fact we would like that for small = AGrrii, the results do not differ much 
from the flat metric, i.e. the metric g^^, is just slightly perturbed from the flat 
metric 77^,^. This fact and the fact that {B~^ez)°' must produce a cut in the z- 
plane across which a vector must rotate over an angle 2nai determines the local 
exponents to be: 

/ 



A? 



-a, I (61) 



See also |17|. Now we must take G (0,1) in order to fuUfill relation 
Also if we would exceed this limit the total mass would become too large and the 
universe would close. 

Note that if we would transform to a different global Lorentzframe, using a 
Lorentztransformation A, the effect on the monodromy would be 

M„ A-^M„A (62) 

This is just a conjugation that doesn't change the structure constants of the group 
(S0(2,l)). So the abstract group remains the same. Finally it is important to 
remark that if we follow a trajectory around all particles in the positive direction 
(anti-clockwise), this is equivalent to follow a trajectory around infinity in the 
negative direction, i.e.: 

Mfot = B^R^B^ = B^R^B'^^ BiRiB^^ (63) 

We will treat infinity as a singular point on the Riemann sphere, just as the par- 
ticle positions. However because we demand that our universe is open the nature 
of this singularity will be slightly different from the particle singularities. We are 
now ready to reformulate our problem of solving the e/ as an old mathematical 
problem, the Riemann-Hilbert problem: 

Given the following transformations for parallel transporting a vector around a 
closed loop: 

e/ - (64) 
[l]n eni{C - {ai, ...ttN, 00}; Zo) 

M^e 30(2,1) 

Find the multivalued, holomorphic (=meromorphic) functions e/ that transform 
in this way, with initial value eJ^{zo) =e^ . 

It is actually more convenient to work in the locally isomorphic group SU(1,1). 
Denote C°'{z) to be a spinor, transforming under SU(1,1) and C*° to be a spinor 
transforming under the complex conjugate representation of SU(1,1). (Note that 
is not the complex conjugate of C°-) If the components of are given by 
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then the spinor transforming under the complex conjugate representation is: 




(66) 

where /3 is an arbitrary real constant. The vector e/ transforming under S0(2,l) 
is now constructed as follows: 

/ 2ab \ 

e'^ = /5 -iV2a^ (67) 
V tV2b^ J 

Note that we automatically construct a null-vector e/ in this way. [| In the fol- 
lowing we will always use this spinor representation SU(1,1) of the Lorentzgroup 
in 2+1 dimensions. The (modified) Riemann-Hilbert problem we have to solve is 
thus: 

C M: (69) 
Hn e 7ri(C - {ai...aN, oo}, Zo) 
M„ G S[/(l,l) 

The Riemann-Hilbert problem has a long and interesting history. Riemann was 
the first (1857) to work with a system of functions satisfying a linear differential 
equation with a given monodromy group (however, he couldn't prove that this 
system of functions exists). In 1912 Hilbert added this problem to his list and 
since then many people worked on the problem (and solved it; that is, proved 
that there is a system of functions with a given monodromy group). People 
who worked on the problem were: Plemelj [|I2|, Birkhoff, Rohrl (who proved the 
statement for an arbitrary Riemann surface) and Lappo-Danilevsky. Another 
important breakthrough was the reduction of the Riemann problem to a system 
of Schlesinger equations. For a review on this see 0, We will follow in this 



paper two important ways to solve the Riemann-Hilbert problem. First we will 
use a second order Fuchsian differential equation to attack the problem and solve 
the simple case of 2 particles. Then we will use the approach of a system of first 
order differential equations and the Schlesinger deformation equations to write 
down a general (formal) solution first found by Miwa, Sato and Jimbo [|1^ . 



^The fact that we can express the components of C* in the components of C is a consequence 
of the fact that for SU(1,1) the conjugate representation is not independent but can be written 
as: 

-1.. / 1 
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3.1 The Fuchsian differential equation 

It was already recognized by Riemann that differential equations with regular 
singularities have a monodromy group. If one would take the d independent 
solutions to the d'th order linear differential equation with N singular points 
(infinity is not included in the number N, but it is singular) and move them 
around a singular point they would transform into a linear combination of them. 
The general form of the differential equation with only regular singular points 
(i.e. Fuchsian singularities, see below) was found by Fuchs and reads: 

gi,)^ + ^(A^-l)(^) g(d-l)^ + +.... = (70) 

Q{z) Q\z) 

P;(7v-i) is a polynomial of order 1{N — 1) and Q = Y[k=i{^ ~ o,k)- 
As we will use the 2 dimensional representation of the monodromy group SU(1,1) 
we can take d=2. If we had chosen to work with SO (2,1) we would have to take 
d=3. We want that the 2 independent solutions and to this equation are 
"regular" singular near the points a^. This implies that by definition they must 
have the following behaviour near ak- 

C {z - akY^\hl + cl{z - a,) + ...] (71) 

This local behaviour is defined in a local basis where the monodromy is diagonal. 
In general the solution near a singular point is a linear combination of the 2 
solutions [a = 1,2) given above (see (^5))). The coeficients (called the local 
exponents) must obey the Fuchs relation: 



EEA^ + EA:. = ^^^(iV-l) (72) 



d N d 
=lfc=l a=l 



The first non-trivial case appears with 2 singularities (and infinity). This is also 
a special case because as we will explicitely show, the local exponents completely 
determine the Fuchsian differential equation and its solutions. With N singular 
points there will be (N-2) undetermined parameters (free parameters) in the 
differential equation. They must be determined by matching the monodromy 
matrices M„ with the general monodromy of the solutions to that equation. If 
we have more than 2 singularities a complication is that apparent singularities 
may occur, (at most N-2 [|l^). These are singular points in the differential 



equation but not in the solution. They will however generate zero's in the At 
a zero the local exponents are possitive integers, and they are counted in the sum 
(|72[). This is reason why for more than 2 particles the method of finding a d'th 
order Fuchsian differential equation becomes very difficult and we use another 
method. Another remark is that the Riemann- Hilbert problem is not uniquely 
solvable if we only give the monodromy matrices. Let us define (f){z) to be the 
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d independent solutions to a definit Fuchsian differental equation. For a given 
monodromy there will exist d such sets of independent functions <f>i{z), which are 
solutions to d Fuchsian differential equations with the same monodromy group 
but with different exponents at the singularities. These d sets are rationally 
independent with respect to eachother. I.e. a new set can always be written as: 

m=^2q^iz)■Mz) (73) 

i=l 

and qi{z) are rational functions, (rational functions will not change the mon- 
odromy of the solutions). A basic result is ||15|| : 

Only if we specify the local exponents exactly (so not only up to an integer 
that would not change the monodromy, but also the integer part) then we can 
find a unique solution to the Riemann-Hilbert problem. As it turns out, the lo- 
cal exponents at are proportional to the masses of the particles and the local 
exponent at infinity to the total mass. As we saw in section 3, physics demands 
that we choose the masses in the range (0, 1). This determines the integer part 
of the exponents and makes the problem uniquely solvable. This method of find- 
ing a Fuchsian differential equation that has the right monodromy is particularly 
suited for the 2 particle case. We will use a different strategy for the general case. 
Below we will work out the 2 particle solution. 

3.2 The 2 particle solution 

To solve the 2 particle case one has to notice that a second order Fuchsian dif- 
ferential equation with 3 regular singularities, one of which is located at infinity, 
can always be written in the following form: 

d'C+j: ^' dC+{j: - "^-^^^ +ALAL}-^ r = (74) 

where we have defined ajt cyclic: as = etc. The general solution can be written 
symbolically using the Riemann-scheme: 

{Oi 02 CXO \ 
\\ \\ -.zS (75) 

A? \l XI J 

From this notation one can find in the upper row the singular points and in the 
second and third row the local exponents. Remember that these local exponents 
have to satisfy the Fuchs-relation: 

E A^ + E A^ = 1 (76) 

k,a=l a=l 
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Because we still had conformal freedom left in our gauge we can bring the points 
{ofe} to respectively {0, 1} using the following conformal transformation: 



w 



{z - ai) 
(02 - ai) 



(77) 



We would like to transform this equation to the hypcrgcomctric equation. This 
is possible by a transformation to a new variable ^{w), defined by: 



The Riemann-scheme now looks like: 

1 00 

- A2 - A2 A^ + A2 + 



(78) 



w 



(79) 



And this is precisely the Riemann-scheme for the hypergeometric function F{a, b, c; w) 
with coeficients: 



a 
b 
c 



Ai + A^ + AL 
Ai + A^ + A^^ 
l + Xl-Xj 



(80) 



There exist 2 independent solutions to the hypergeometric equation near w 
(the series converges in the circle l^l < 1): 

= F{a,b,c;w) 



w 



1-c 



F(l + a - c, 6 - c + 1, 2 - c; w) 



(81) 
(82) 



And also 2 rationally independent ones which one obtains by taking A^ = A^ + 1 

and A^ = A^ — 1. We will however use the first ones, because they provide us 
with the right asymptotic behaviour. This means that the solution for around 
w — Q becomes (using some basic properties of F): 



C = w^'^{l-w)^^^F{l-a,l-b,2-c;w) 



(83) 
(84) 



We will call this solution = (P(Aj), P(A^)). The argument of P determines 
the local exponent of that solution. If we want to find the solutions outside the 
circle \w\ — 1 we have to take the analytic continuation of this solution to the 
rest of the complex plane. Of course this is all well known for the hypergeometric 
equation. Denoting the 2 independent solutions to the Papperitz equations near 
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w = 1 and = oo in the way as above we find that the analytic continuation 
looks like: 

P{Xl) = AlP{Xl) + AlP{Xl) 
P{Xl) = BlP{Xl) + BlP{Xl) 
P{>^\) = AlPiXl) + AlPiXl) 

^A?) = BlP{Xl) + BlP{Xl) (85) 



The coefficients can for instance be found in |16]. Using the local expansion ([71|) 
we can simply see what happens if we move around the location = {0, 1, oo}. 
The -P(A^) will change as follows: 



PiXt) - e'-^'^^PiXt) (86) 



we can actually derive the monodromy explicitely ([0). In the following we will 
give these monodromy matrices and match them with the monodromy we impose 
for a particle. First we define the following basis: 

u = RAP{Xl) (87) 
V = S P{Xl) (88) 

Here R and S are functions of the parameters of the system only and A is: 

r(i + A? - x\nx\ + XI + xinx\ + XI + XI 



A 



r(i + A} - Af)r(Af + Ai + A^)r(Af + Ai + a^ 

r(i + xj- x\)r{x\ + XI + xl)r{x\ + xl + x 



oo/ 

' ) 

oo/ 



r(i + A} - A?)r(A? + XI + xmxi + XI + xi) 

Now the monodromy for moving around = is: 



(89) 





And the monodromy for moving around w = 1 is: 



^1 = n „2.a? (90) 




/ 



Mo 



S (e 2-e 2) fxe 2 - 



(91) 



^ sin^(Al + A^ + XI) sin7r(Af + A^ + Xl) 

^ sin4A2 + Ai + Aysinvr(Al + Ai + Ay ^ ^ 

This monodromy has to be matched with the desired SU(1,1) monodromy. We 
will fix particle 1 (at u; = 0) and allow particle 2 (at = 1) to move in an 
arbitrary direction with an arbitrary velocity (see Fig. 3). The monodromy for 
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• • • 

z=0 z=l z=oo 

Figure 3: Two particles moving on the complex plane 
particle 1 is then: 

Ri=i e-"i ) = ^^""^ ^^^^ 

The second particle has a more general monodromy of the form: 

^ ^ _ / cos(7ra2) — ^ sin(7ra2) coshr] — sin(7rQ;2) sinh ?]e*^ \ 
^ ^ ^ I — sin(7rQ;2) sinhr^e"*^ cos(7ra2) + ^sin(7rQ;2) cosh?] j 

(94) 

The monodromy around = cxd is now fixed by the relation: 

M1M2M00 = / (95) 

Matching the monodromies gives: 

2X1 — (mod rii) 

2X1 = '^1 (mod 712) 
2A2 = — tt2 (mod 713) 

2A2 = a2 (mod 71,4) (96) 
= coth2(i?7) (97) 
i? i 

- = -sinhr^e*'' (98) 
o Z 

We take here (mod Ui) rii G Z, instead of (mod 2ni). This has to do with the 
fact that SU(1,1) is the covering group of S0(2,l). In this case this implies both 
the transformations ±5" e SU(1,1) correspond to one element in S0(2,l). For 
reasons explained before (section 3) we will take = (See also ^^). We can 
simplify relation ( |97D to: 

cos(27rA^) = cos(27rA}) cos(27rA2) - cosh r] sin(27rAi) sin(27rA2) (99) 

This is the precisely the same relation found in p[. We will therefore interpret 
47rA^ as the total deficit angle of the universe which means that atot = 2A^ must 
be in the range (0,1). Next we can use the Fuchs relation ( [76D to find A^: 

Xl = l-Xl (100) 
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The complex conjugate representation can be retrieved by using (p^). Now we 
can write down the full solution to the monodromy problem: 



sinh ?7e*''AP(— «!, 0:2, atot', w) 



y P{ai,a2,atot;w) 



(101) 



with: 



P{ai,a2,atot;w) = w ^ {1-w) ^ F(l + y + ^ ^, y + y + — , 1 + ai; «;) 



Using the formula (|67D we can write down the dreibein: 



/3 



/ ismhr]e'^AP{ai)P{-ai) 
^sinh^ r]e^'^A^P\-ai) 
[ tV2P\a,) 



(102) 



(103) 



It is now interesting to check the asymptotic behaviour found in formula (|101|) . 
Since we know the local exponents at infinity we can write down the following 
expansion for P(A^): 



w 



w 



+ 



with 



atot 



(104) 

(105) 
(106) 



Using (^TD we find the asymptotic behaviour for ej^: 



I 



w~ 
w~ 



-atot 



\ W 



atot-'i 



W 



+ 



(107) 



It is important to realize that the global Lorentzframe we use to calculate the 
asymptotic behaviour is different from the one we used to calculate the P(A"). 
In this frame particle 1 is not moving and as a result Mi is diagonal in the basis 
P{Xi). In the frame we used to derive the asymptotic behaviour the monodromy 
around infinity is diagonal which means that the effective center of mass is not 
moving. The relation between the 2 frames is given by the coefficients and 
of ([85|). We see that the asymptotic behaviour is the same as in ( |i6|Ji7|P8| ) . 
One could also try to find the total angular momentum by matching the first 
coefficient in the expansion with ( |46|J47|P8| ) : 



Bt. 



(108) 
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3.3 the N-body case 

To treat the N-body case we will use a completely different approach, first followed 
by Schlesinger. Every monodromy problem will correspond to a first order matrix 

differential equation. First we write C"(-2) = '^'^isC^i^o) where Y{z) will be a 
fundamental matrix of solutions to the differential equation. The Y[z) has the 
following properties: 

i) Y{zo) = / 

ii) Y{z) is holomorphic in (C — {ai...a]y, oo}) 

iii) Y{z) has the following short distance behaviour near a singular point: 

y(^)^~"'(^-a,)^'y,(^) (109) 

where Mj = e^'^*^' and Mj are the monodromy matrices. It is also important that 
the matrix Yi{z) is holomorphic and invertable at 2; = Oj. So: 

Yi{z) (y°(a,) + Y\a,){z - a,) + ....) (110) 

with det|y°(ai)| ^ 0. 

To derive the form of the differential equation we first consider the object Y~^{dY). 
This is also a single valued expression (the monodromy cancels out) and holo- 
morphic at (C — {ai...aN, 00}). By using the local behaviour near the singularity 
Oj we can find an expression near the location af. 

Y-\dY) ^-"^ Yr\z)^^Yi{z) + holomorphic at z = (111) 



Hence it must be of the form: 



Y-\dY) = J2 (112) 

with 

A, = Res,^ay~\dY) = {Y^Y^UY^ (113) 

The behaviour of this equation at infinity can be studied by transforming to a 
new coordinate w = -. We find that that: 

z 

N 



Aoo=-EA (114) 



To prove that the solution to the differential equation: 



N 4. 

dY = Y,y—^ (115) 
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is unique we consider 2 solutions Yi and Y2 and study the function YiY2^^. Using 
the conditions ii) and iii) we see that this function is single valued and holomor- 
phic everywhere on the Riemann sphere. This implies that it must be constant. 
Using i) we find that it must be actually / and we have proved the statement: 
Yi = Y2, which in turn proves that the solution is unique. The proof that there 
actually exists a solution is much more involved and was only esthablished for 
the first time in its full generality by Rohrl. We will state the result obtained 
by Lappo-Danilevsky. He obtained a formal solution in terms of expansions in 
hyperlogarithms. See for a review. He proved that if the Lj were sufficiently 
small to ensure the convergence of this series, then the Riemann problem has 
d rationally independent solutions. Every other solution is given as a rational 
combination of these solutions. Remember that if we multiply Y with a rational 
matrix R{z) the monodromy does not change. Even when no convergence was 
guaranteed the solution could still be interpreted as analytic continuations. If 
we furthermore specify the local exponents exactly, that is specify the local be- 
haviour of the spinor exactly (in a canonical basis), then we can find a unique 
solution for up to an overall constant. This constant must then be determined 
in our case by matching the series expansion of e/ at infinity with (^7|) . An 
important additional condition is that the eigenvalues of Ai, An, A^o (which 
are equal to the eigenvalues of Li, L^, Loo) do not differ by integers. In that 
case logarithmic singularities appear which we will not include in this paper. It 



was Miwa, Sato and Jimbo fl^ who wrote down a solution to this monodromy 



problem in terms of correlators of free fermions with twistoperators at the posi- 



tions of the singularities. It was Moore [IS] who translated it in the language of 



a conformal field theory for free fermions. Their solution was: 

^ {Voo{oo)VN{aN)...Vi{ai)) 

The operators Vi are twistoperators to be defined shortly ( |126| ). First of all, it 
is easily seen that taking the limit 2; ^ we can write 



(Voo (00) Vjv (ajv ) . . . Vi (ai ) ) 

-|- holomorphic terms ai z = Zo 

Where we have denoted the Wick-contraction as: 



r{z)i^p{zo) = — ^ (118) 



^Note that this relation aUows in principle for a more general set of fields and cp{zo) 

which have spin A and 1 — A respectively (a b-c system). We want however that the spin 
(conformal weight) of C"(^) is the same as that of &"(z), so that they transform the same way 
under local transformations of the coordinate z. This gives that A = ^ as we have chosen in 
the text. 
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Figure 4: N Particles moving on the Riemann-sphere 
Using the above formula we find that: 

yp{zo) = 5^ (119) 



as required by condition i). An important difference with [|18l is that we have 
put a twistoperator at 2; = oo, altering the "outgoing vacuum" and changing 
the scaling dimensions of the fields at finite distances Oj. The twist at infinity 
must reproduce the monodromy at infinity and contains the information about 
the total mass of the universe. 

We now come to the subject of choosing the twistoperators. They must be 
constructed in such a way that they reproduce the correct monodromy properties 
at z = ttj. In the following we will consider a particle at z = with boost-matrix 
Bi = B and twistoperator Vi(0) = V{Q). We will drop the index i in the following. 
The same story holds for all the particles. For the particle at 2; = we want: 

^"(e^-*^) = M^'p^^^z) = {BRB-^Ypij''{z) (120) 

This means that the short distance behaviour of ifj'^lz) and V^(0) must be of the 
form: 

ViQ)riz) = Y^iB^Tf^z^'^H^iz) (121) 

where B is a boost-matrix, Dp are the eigenvalues of the matrix L and is a 
holomorphic vector at z = 0. We can write this as: 

J2 B'^z^'^H^ = J2iB^iz)TpH^ (122) 

/3 /3 



where we supressed the index i and where: 

rt 
2t 



n{z)=l^' ° ) a = 4Gm (123) 
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Due to Q{z) this expression is multi-valued and, by analytic continuation, we find 
if we move around the particle at z = af 

{BVL{z)YpH'^ {BRB-^Y^{BVl{z)HY (124) 

where: 

/ -i-wa Q \ 

«=( A) (125) 

This is of course precisely the right monodromy and our task is to find the 
twistoprators with this property. In ref.[18| one finds: 

V =: exp[-p, / dw i>aCp^pp] : (126) 
J c 

where : ilJa{w)t'^p'i/jfj{w) := J" is the su(l,l) Kac-Moody current, and C is a 
contour along the cut in the plane from the position z = to oo, and —Patap = 
Laf3. In order to obtain the operator product expansion (OPE) at 2; = we 
note that we can use the freedom to transform the Dirac-spinors -0" by a global 
SU(1,1) transformation. 

^°(z) ^ 7/;° = A^'f^^^iz) (127) 

Next we write: 

V{0)^P''{z) = : exp[ / dw VL^^uV] ■ ^°(^) (128) 



c 



Defining: 



exp[/" dw {il)BY{B-^LB)^YB'^i^y] : ^/'"(z) (129) 
J c 

V!" = (5-V)" (130) 

I" = mr (131) 



we find: 



V{QW{z) =■ exp[ / dw : 5"^^^(^) (132) 

J c 

Here D^^, = Pot^^^ = fi°,^, « = 4Gm and t°^,=diag(l,-l).The next step is to 

bosonize the -0", ^jJ and V^(0) (= V"(0)). Of course we can not simultaneously 
diagonalize all the Li, but we don't need that actually because we only need the 
OPE locally at 2; = 0. This result is independent of the representation that we 
choose (i.e. fermionic or bosonic). So we write: 

^"(z) = : e*^"(") : (133) 
^(z) = : e-^^°(^) : (134) 
Y^-.^'^dJ'^: = Y.DadVa = D-d^ (135) 

a a 

y(0) = : e'^-^(°) : (136) 
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Where Da are the eigenvalues of L. Formally we should add Klein-factors (c^ and 
Co) to the expression : e**^" : and : e"**^" : to ensure that the fermions anticommute 
for all a. We will leave them out in the following because we won't need them 
futher. Note that if we bosonize the twistoperators ( |126| ) according to the rule 
([135D , we also pick up a term at infinity: 1/(0) = e'^''^^^^ e~'^-'^^'^\ This extra 
term is of course part of the twistoperator at infinity. There is thus a slight 
inconsistency in the notation. Because the twistoperators are really nonlocal 
operators they can not be assigned to 1 point. 
Using the bosonic representation and the following identity: 

. ^it^-viz) ■■ Qi>'-V{w) ._ _ ^^K-^ ■ QiK-v{z)+i\-0{w) . (137) 

we find for the OPE: 

\/(0)^"(2) = J2 B'^^z^^H^ (138) 

where is again a holomorphic vector at 2; = 0. This is precisely expression 
(|121|) . So we esthablished the right behaviour of the solutions near the points 



Oj. The next question is obviously to check the behaviour of this solution at 
infinity. To find that one needs to know how the field ip°'{z) scales, i.e. we need 
to know the scaling dimensions A"" = 2/i" of the fermion. However the fact that 
we have put a twist operator at infinity spoils the naive scaling dimensions (see 
[0). Before we do that we like to get rid of the Zo dependence. More precisely 
we want to replace condition i) in the beginning of this section by a different 
condition at 2; = 00. That new condition is: 

i') 

= I (139) 

To implement this new condition we first have to bring the field i^f^izo) to infinity, 
using the rule: 

^^3(00) = lim M^o)z^^ Af, = 1 (140) 



In front of expression ( 116 ) we have a factor {zo — z) and we use this factor 



{{zo — z) Zo if Zo 00) to take the limit ( |140| ) . However there is already a 
twistoperator at infinity which implies that we will have to deal with another 
divergent factor. For convenience we will assume that the center of mass is not 
moving, which means that at infinity the monodromy is diagonal. The twist 
at infinity is defined by the total action of the twists at finite points. Because 
all these twists contain integrals from the particle position ai to infinity, and at 
infinity all these twists do not add up to zero, there is a nontrivial twist at infinity, 
see Fig. 4. 

To see what happens if we bring the operator ipf^lzo) close to this twistoperator 
we first bring Voo to a point b in the finite plane and bosonize: 

V^{b)Mzo) =: e-'"^--^^'^ :: e'^^^^^") :~ {b - Zof^Hp{z) (141) 
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where Hf^^z) is holomorphic at z = b. This suggests that if we want to take the 
hmit Zo ^ b we have to do the following: 



lim (6 - V^{b)Mzo) =: V^{b)ij^{b) : (142) 

Zo-*b 



One can check that if we take the limit: 

-D 



\im{b- Zo^^YI (143) 

we have precisely the condition i') at z=b: 

Yb = I (144) 

This means that in taking the fermion ipi3{zo) to infinity we need to multiply by 
a factor Zo ^ in order to obey the new condition i') at infinity. All together we 
write now: 



{VooVn...Vi) 



To determine the asymptotic behaviour of Y' we have to calculate the scaling 
behaviour of the iIj°'{z). To do that we again bosonize and note once more that 
we have chosen the center of mass to be at rest, meaning that \4o has a "diagonal 
form". As mentioned before, the operator : K^o'^/j : at infinity spoils the naive 
scaling behaviour. In the bosonized picture we have the following operator at 
infinity: 

Of =: e-*^--^-^v'/3 : (146) 



Following (||T9|) we must change the scaling dimensions by the following rule: 

A'" = 1 - {D^ + 5^) (147) 
This means that the asymptotic behaviour of Y' looks like: 

(148) 



z 2 nu[z) z 2 -n^2 

Z~^''^H2liz) Z^H22(Z 



where Hij{z) are holomorphic functions as z ^ oo. Next we choose C = (1, 0) so 
that we have: 



- "tot 



C'-^{ '^^1 + + (149) 



1 



Z 2 



This is in fact the right asymptotic behaviour at ^ — oo. One can actually 
find an explicit expression for these correlators in ([|T5). This is a rather formal 
infinite series and hardly of any practical use to us. Therefore we will not repeat 
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it here. The first order matrix differential equation can be viewed as an equation 
for parallel transport on a Riemann-sphere. We can define a connection: 

A = Y.^ (150) 

and a covariant derivative 

V=^^-A (151) 

In this picture e/ is determined by parallel transporting the e"{zo) from Zo to 
z. The result depends on the path we choose. Because the connection A is 
holomorphic everywhere except at the singularities and oo the result will only 
depend on the first homotopy class of the punctured sphere. In general the 
dependence of the Ai on the points can be very difficult. If we use normalisation 
condition i) then also a dependence on Zo is possible. Schlesinger determined 
a differential equation for this dependence by demanding that if we move the 
particles to another location (for instance by time-evolution!) the monodromy 
will not change. This is because we don't allow the particles to cross a string 
connected to one of the other particles, by smoothly deforming the string to 
another position as a particle comes close to a string. This also means that 
the topology of the strings will become very important in the dynamics of the 
particles. Another way of putting the fact that the monodromy does not change 
during time-evolution is that the fiat momenta pa are conserved. They act as 
a kind of charges in this picture. The equations found by Schlesinger are called 
isomonodromic deformation equations and are given with the normalisation i') 
by: a 

-r (152) 



dai z — ai 

dA ^ [A„A,] 



i^j (153) 



dA ^ _^[AM (^54) 



In order to get a feeling for the solution presented in ( |116| ) we will calculate a 



simple example, namely the case where all Ai commute. From ( |152| , |153|J154|) we 



see that the Ai do not depend on the particle's positions a^. Next we diagonalise 
all Ai simultaneously. This means that we are describing a static configuration 
of particles. It is very easy to find the exact solution to the differential equation 



y = X{i^r (155) 



_ -j^ Zq Oil 



^If we choose normalisation i), the equations will also include a Zq dependence, see |Q and 



27 



where : 

= ( T ^ ) = ^^""^ ^^^^^ 

The answer in the normahsation i') is according to ( |145| ): 

N 

Y' = lirn^ z^°^Y = Y[{z- a,i)^^ (157) 

° i=l 

In the static case it is indeed true that Mtot = Y^iLi ^i- Next we choose the 
vector C = (1, 0) so that we find: 

N 

r = (n(^-a.)-^,0) (158) 

i=l 

and 

N 

e/ = (0,n(^-«^)~"%0) (159) 

i=l 

Now using the formula (^9]) we can calculate the conformal factor defined by 

= e1{z)e,^a{z) (see (P): 

TV 

= X{{z-ai){z-ai)-^^ (160) 

which is the correct answer (see 0) On the other hand we would like to calculate 
the correlator ( |145| ). Because all Li commute we can bosonize the fermions at all 
points Oj simulteneously. We may view the fermion 'il)"'[z) as two twistoperators 
with "charge" 1 and ^'/^(^o) as 2 twistoperators with "charge" —1. The correlators 
that we have to calculate (a, = 1, 2) are only nonzero if the total charge adds up 
to atot- This is due to the infrared cut-off that has to be taken care off |T^. Only 
when the total charge is atot we find a nonzero correlator. Now using repeatedly 
formula (|137|) and noticing that only ioi a = j3 the total charge adds up to atot 
we find: 

r^ite^ii^fe^^ni-'..^ ,161) 

\.\.i<j\^i ~ (^j) ^ i=l 

The use of the operators at infinity is now to ensure that the total charge adds 
up to atot- 

N 

^ai = atot a = (3 (162) 

i=l 

We could call this a Coulomb-gas picture since it bears much resemblance with 
the background Coulomb-gas of ref.lH]. Finally we want to say something about 
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the correlator: {VoqVn ■ ■ -Vi) . We can derive a kind of Knizhnik-Zamalodchikov 
equation for this correlator (see also [l^). For that consider first: 



(163) 



{VooVn...Vi) 

) 

(164) 



S'^ {V^f^VN...Vi) 



where we have defined: 

1 - d I d - 

Tpa =■ ^i^p{Zo)^^a{Zo) " l^-Q^i^ p{Zo)tpa{Zo) ■ (165) 

Because we have a background charge at infinity, the naive "energy- momentum 
tensor" T = trT has to be corrected by a boundary term: 

T = f-Y. \D^d : T^a^a : (166) 

The boundary term however does not change the simple pole term in the OPE 
so that we still have: 

=um=i (167) 

Jai 2m Jai 27rZ Otti 

So we find: 

iy^fvN...Vi) 1 d 



{V^Vm...Vi) (168) 



/a, 2m{V^VN..yi) {V^Vn...Vi) da^ 
On the other hand we may also write the correlator as: 



+ t:- 



{VooVn---Vi) 2dzZo — z 2dzoZo — z 

(169) 

This can be worked out by using an expansion around Zo'. 

Y = I + A{zo){z - Zo) + ll^A^o) + A\zo)]{z - Zof (170) 

A is defined in ( [1501 ) . Integrating this expression around z^, = ai and comparing 
with ( p.6(j|) we get: 

^iy^v^...v,) = \i dzo E , f;^^' , (171) 



•"in the mathematical hterature this is caUed the tau-function. 
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As mentioned before there also exists a Schlesinger equation in which Zo is in- 
cluded in the set deformation parameters (see fl^)- The equation we need is: 

^ ^ [^ii ^j] {112) 

OZq j Zq Qjj 

We can use this formula to find that tr AiAj does not depend on Zo at all! This 



allows us to calculate the integral in ( |171| ) and we finally find: 



^{V^V^...V,) = lT.'i^{yooVN...V.) (173) 

Although this looks a bit like a Knizhnik-Zamalodchikov equation, there is are 
differences. Most importantly, the term tr A^Aj may depend in a very compli- 
cated way on the positions of the particles Oj. The fact that we didn't find a 
Knizhnik-Zamalodchikov-equation (with constant matrices Ai) is interpreted in 
as a result of the fact that on the Cauchy-surface the translational symmetry 
is broken. 

As solving the Riemann-Hilbert problem is equivalent to finding the Ai, we have 
in general not an explicit expression for them. Only in some special cases, where 
we are able to find the Ai we can write down a full solution to the Riemann-Hilbert 
problem. It is also very interesting to look at the OPE of 2 twistoperators and 
determine the monodromy if one particle moves around the other. Alternatively 
one can look at the exchange-algebra (or half-mono dromy) and determine what 
happens if 2 particles are exchanged. We expect that these half-mono dromies 
obey some non-abelian representation of the Yang-Baxter equation (see 
We hope to say more about this in the future. 



4 Discussion 

In a series of papers @, 0] 't Hooft solved the classical N-particle problem in 
2+1 dimensional gravity. In his approach he used flat coordinates that obey 
unusual boundary conditions. As a result these coordinates become multivalued. 
In this polygon approach he developed a Hamiltonian formulation for the phase 
space variables^ and by replacing Poisson-Brackets by commutators he was able 
to formulate a first quantized theory. It is however rather difficult to write down 
an explicit wave function in this formalism as all kinds of transitions that could 
take place have to be taken into account as boundary conditions on the wave 
functions. Surely if one wants to formulate a second quantized theory these 
boundary conditions become untractable. This was the main motivation for us 
to define a single- valued coordinate system, much like in the case of anyons where 

^-^The phase space variables are the lengths of the edges of the polygons and the rapidities 
(boost-parameters) defined across the edges. 
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one has multi-valued and single-valued gauges. In the case of anyons the single- 
valuedness has a price, namely the introduction of an electromagnetic interaction 
between the anyons. Also in the case of 2-1-1 dimensional gravity we can trade 
the multivaluedness for a nontrivial gravitational interaction. This implies that 
we have to introduce a nontrivial metric g^u 7^ Vtiv 

The physical quantities of the theory are the particle positions and the con- 
jugate momenta. As the gravitational field doesn't have any degrees of freedom 
it is in principle possible to integrate them out. To do that one must first study 
the gravitational field in our particular gauge. That is precisely what we did in 
this article and we found that the 2 particle case could be solved explicitly while 
in the N-particle case, the gravitational field (or rather the dreibeins) could be 
represented in terms of correlators of a two dimensional conformal field theory. 

In order to reduce the system to its physical degrees of freedom we need to 
find a Hamiltonian. In section 1 it was argued that the Hamiltonian is given by 
a boundary term: 

H = J d^x Vh^^^R= <f dz d(f){zz) + dz d(p{z,z). (174) 

where is defined by: e'^ = e1{z)ez,a{z). This expression can be simply evaluated 
by using the monodromy-matrix at infinity and the asymptotic behaviour of 
(^2|). One finds the expected result; H = atot, the total deficit angle of 



[z,z 



the universe. By using (p5|) it then easy to find the well known result: 



N 



27r2i/t°^ = log[ne'"^1„/3 (175) 



i=l 



where Li = —p\t1^p ; p^a^ are the flat momenta and are the generators of 
S0(2,l) (or equivalently SU(1,1)). Because we restricted ourselves in this paper 
to particles with positive mass, and didn't consider anti-particles with negative 
mass (which could in principle be defined) we have positive total energy. It is 
however also particular simple in our gauge to prove that for arbitrary matter 
couplings (that obey reasonable boundary conditions at infinity) the total energy 
of an open universe is positive definit (see 0, |11[). This is of course a nessecary 



property in order to have stable vacuum in a second quantized theory. It is 
also important to realise that gravity in 2-1-1 dimensions contains long range 
interactions. This can be seen most easily in the fiat coordinate-system because 
the wedge cut out of space-time extends out to infinity. This means that we don't 
have an asymptotic fiat region. If we want to define an S-matrix in the future we 
know that we will have to define in- and out-states with great care. We hope to 
say more about the particle dynamics and the quantisation of this system in the 
future. 
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Note added 



While writing this paper I found out that A. Bellini et al. [1^ published a preprint 
concerning the same issues we adress. 
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